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Charge carrier transport in single-layer graphene with one-dimensional charged defects is studied 
theoretically. Extended charged defects, considered an important factor for mobility degradation 
in chemically-vapor-deposited graphene, are described by a self-consistent Thomas-Fermi potential. 
A numerical study of electronic transport is performed by means of a time-dependent real-space 
Kubo approach in honeycomb lattices containing millions of carbon atoms, capturing the linear 
response of realistic size systems in the highly disordered regime. Our numerical calculations are 
complemented with a kinetic transport theory describing charge transport in the weak scattering 
limit. The semiclassical transport lifetimes are obtained by computing scattered amplitudes within 
the second Born approximation. The transport electron-hole asymmetry found in the semiclassical 
approach is consistent with the Kubo calculations. In the strong scattering regime, the conductivity 
is found to be a sublinear function of electronic density and weakly dependent on the Thomas-Fermi 
screening wavelength. We attribute this atypical behavior to the extended nature of one-dimensional 
charged defects. Our results are consistent with recent experimental reports. 

PACS numbers: 81.05.ue, 72.80.Vp, 72.10.Fk 
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I. INTRODUCTION 



The isolation of graphene — the queen of two- 
dimensional materials due to its remarkable physical 
properties — by the exfoliation method has triggered in- 
tensive studies of its fundamental properties and has 
opened horizons for future technologies. 1_3 Since that 
time, various methods of graphene growth have been ex- 
plored in order to make the fabrication process scalable; 
a prerequisite for developing graphene-based devices and 
technologies i 4 Nowadays, several techniques are capable 
of producing high-quality, large-scale graphene. These 
include epitaxial graphene growth on SiC, 5 and chemical 
vapor deposition (CVD) of graphene on transition metal 
surfaces^- The advantages of the latter method lie in its 
low cost, possibility to grow large graphene sheets (tens 
of inches), and ease of its transfer into other substrates »£ 
Currently, there is a strong motivation for exploring elec- 
tronic and transport properties of CVD grown graphene 
because it represents one of the most promising materials 
for flexible and transparent electronics. 

The studies of the transport properties of graphene 
are often focused on the fundamental question: what 
limits a charge carrier mobility in it? As far as CVD- 
grown graphene is concerned, it is believed that its trans- 
port properties are strongly affected by the presence of 
charged line defects^ Usually, the growth of graphene 



by the CVD-mcthod requires to use metal surfaces with 
hexagonal symmetry, such as the (111) surface of cubic 
or the (0001) surface of hexagonal crystals^ The mis- 
match between the metal-substrate and graphene causes 
the strains in the latter, reconstructs the chemical bonds 
between the carbon atoms and results in formation of 
two-dimensional (2D) domains of different crystal ori- 
entations separated by one-dimensional defects^— The 
nucleation of the graphene phase takes place simulta- 
neously at different places, which leads to the forma- 
tion of independent 2D domains matching corresponding 
grains in the substrate. A line defect appears when two 
graphene grains with different orientations coalesce; the 
stronger the interaction between graphene and the sub- 
strate, the more energetically preferable the formation of 
line defects is. These line defects accommodate localized 
states trapping the electrons, originating lines of immo- 
bile charges that scatter the Dirac fermions in graphene. 

It is well established that the presence of grains and 
grain boundaries in three-dimensional polycrystallinc 
materials can strongly affect their electronic and trans- 
port properties. Hence, in principle, the role of such 
structures in 2D materials, such as graphene, can be even 
more important because even a single line defect can 
divide and disrupt the crystal^ A series of recent con- 
trol experiment o 13 ! 14 strongly indicate that line defects 
are responsible for lower carrier mobility in CVD-grown 
graphene in comparison to the exfoliated samples* 1 ^ - — 



We note in passing that one-dimensional (ID) defects 
have been observed not only in experimental studies on 
CVD growth of graphene films, for instance, on Cu^ 4 - 
NifiS Ir^iS but also in single graphene layer after electron 
irradiation^ and in highly oriented pyrolytic graphite 
surface. 21 Possible applications include: valley filter- 
ing based on scattering off line defects, 22 ferromagnetic 
ordering in line defects ; 23 ! 24 enhancement of electron 
transport^ or chemical reactivity2& due to induced ex- 
tra conducting channels and localized states along the 
line, quantum channels controlled by tuning of the gate 
voltage embedded below the line defect^ and correlated 
magnetic states in extended defects *22, 

Several theoretical studies have been recently reported 
addressing transport properties of graphene with a sin- 
gle graphene boundary 29-31 or polycrystalline graphene 
with many domain boundaries.— On the other hand, 
much less attention has been paid to the effect of charge 
accumulation at these boundaries due to self-doping. 
Transport properties of graphene with ID charged de- 
fects has been studied in Ref. 8 using the Boltzmann 
approach within the first Born approximation. It has 
been demonstrated previously that such approximation 
is not always applicable for the description of electron 
transport in graphene even at finite (non-zero) elec- 
tronic densities i^ 3 - - — In the present work we investigate 
the impact of extended charged defects in the transport 
properties of graphene by an exact numerical approach 
based on the time-dependent real-space quantum Kubo 
metho d 34 i 36 ~ 46 i 48 ' 49 which is especially suited to treat 
large graphene systems with dimensions approaching re- 
alistic systems containing millions of atoms. Our numer- 
ical calculations are complemented with a semi-classical 
treatment going beyond the first Born approximation, de- 
scribing the transport properties in the weak scattering 
regime. 

The paper is organized as follows. The numerical mod- 
els (tight-binding approximation and Kubo approach) 
and obtained results are presented in Sec. |TTJ In Sec. IIIII 
we study the impact of extended charged defects within 
kinetic transport theory. Here, the general expres- 
sion for the scattering amplitude for massless fermions 
within the second Born approximation is derived and 
used to obtain the semiclassical conductivity and the 
transport electron-hole asymmetry. The approaches in 
Sees. HTland lllll provide information about transport dom- 
inated by ID charged defects in distinct regimes. Sec- 
tion IIVI presents the conclusions of our work. Details of 
numerical calculations and analytic derivations are given 
in the Appendixes. 



approach and also present numerical results obtained 
within the framework of these models. 



A. Basics 

To model electron dynamics in graphene, we use a 
standard p-orbital nearest neighbor tight-binding Hamil- 
tonian defined on a honeycomb lattice 3 ^ 50 ' 51 



H = -U ^2 C l C i' + ^ ViC;d 
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where c] and Ci are the standard creation and annihila- 
tion operators acting on a quasiparticle on the site i. The 
summation over i runs over the entire graphene lattice, 
while i' is restricted to the sites next to i; u = 2.7 eV is 
the hopping integral for the neighboring C atoms i and i' 
with distance a = 0.142 nm between them, and Vi is the 
on-site potential describing impurity (defect) scattering. 
Since line defects can be thought as lines of recon- 
structed point defects^r— we model a ID defect as point 
defects oriented along a fixed direction (corresponding 
to the line direction) in the honeycomb lattice. The 
electronic effective potential for a charged line within 
the Thomas-Fermi approximation was first obtained in 
Ref. [fl (see also Appendix A); if there are Mines such 
charged lines in a graphene lattice, the effective scatter- 
ing potential reads as 
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Vj=/ J U-j[-COs(qTFX l j)Ci(qTFX i 

.7 = 1 

^(qTFXij)(Tr/2 - Si(q T FXij))}, (2) 

where Uj is a potential height, Xij is a distance between 
the site i and the j'-th line, <7tf = e 2 kp / '(Treoe r hvp) is 
the Thomas-Fermi wavevector defined by the electron 
Fermi velocity vp = 3ua/(2h) and the Fermi momen- 
tum kp = y/ir\n e \ (related to the electronic carrier den- 
sity n e controlled applying the back-gate voltage). Here, 
— e < denotes the electron charge. The Thomas-Fermi 
wavevector is also commonly expressed as a function of 
graphene's structure constant a g = e 2 kp / (4ireoe r hvp) 
according to <?tf = 4a g A:F- We consider two cases: sym- 
metric, V J; 0, and asymmetric, V > 0, potentials, where 
Uj are chosen randomly in the ranges [—A, A] and [0, A], 
respectively, with A being the maximal potential height. 
In order to simplify numerical calculations, we fit the 
potential ([2]) by the Lorentzian function 
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II. TIGHT-BINDING MODEL AND 

TIME-DEPENDENT REAL-SPACE 

KUBO-GREENWOOD FORMALISM 

In this section, we introduce to the basis of the tight- 
binding approximation as well as the Kubo-Greenwood 



V X =J2 U,(A/(B + CxD) 



(3) 



3=1 



as described in Appendix A. A typical shape of the effec- 
tive potential for both symmetric and asymmetric cases 
is illustrated in Fig. [TJ 
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FIG. 1: (Color online) Effective symmetric [(a),(b)] and asym- 
metric [(c), (d)] potentials describing one of the possible con- 
figurations of 10 (left) and 50 (right) line defects in graphene 
sheet of the size m x n — 1700 x 1000 sites corresponding to 
210 x 210 nm. Maximal potential height A = 0.25it. 



B. Time-dependent real-space Kubo method 

To calculate numerically the dc conductivity a of 
graphene sheets with ID charged defects, the real-space 
order-TV numerical implementation within the Kubo- 
Greenwood formalism is employed, where a is extracted 
from the temporal dynamics of a wave packet governed 
by the time-dependent Schrodinger equation i 32 ' 37-46 ' 48 ! 49 
This is a computationally efficient method scaling with 
a number of atoms in the system TV, and thus allowing 
treating very large graphene sheets containing many mil- 
lions of C atoms. 

A central quantity in the Kubo-Greenwood approach 
is the mean quadratic spreading of the wave packet along 
the rr-direction at the energy E, AX 2 (E,t) = ({X(t) — 
X(0)) 2 ), where X(t) = U\t)XU(t) is the position opera- 
tor in the Heisenberg representation, and U(t) = e~ lHt / h 
is the time-evolution operator. Starting from the Kubo- 
Greenwood formula for the dc conductivity^ 



1-nhp 1 

<r = —^-Tr[v,J(E - H)v x 6(E - H)}, (4) 

where v x is the x-component of the velocity operator, E 
is the Fermi energy, D, is the area of the graphene sheet, 
and factor 2 accounts for the spin degeneracy, the con- 
ductivity can then be expressed as the Einstein relation, 



= e 2 p{E) lim D(E,t), 

t— »OC 



(5) 



where p{E) = pjVt = Tr[d(E - H)]/Sl is the den- 
sity of sates (DOS) per unit area (per spin), and the 
time-dependent diffusion coefficient D(E,t) relates to 



D(E,t) 



(AX 2 (E,t)) 

t 
lTr[(X H (t)-X(0)) 2 5(E-H)\ 



t 



Tr[S(E - H)} 



(6) 



It should be noted that in the present study we 
are interested in the diffusive transport regime when 
the diffusion coefficient reaches its maximum. There- 
fore, following Refs. |4y and |43, we replace in Eq. ([5]) 
limt^oo D(E,t) — > D max (E), such that the dc conduc- 
tivity is defined as 



a = e/p(E)D amx (E). 



(7) 



Note that in most experiments, the conductivity is mea- 
sured as a function of electron density n e . We calculate 
the electron density as n e (E) = n e = J_ oo p(E)dE — 
"ions, where n ions = 3.9 • 10 15 cm -2 is the density of the 
positive ions in the graphene lattice compensating the 
negative charge of the p-electrons [note that for the ideal 
graphene lattice at the neutrality point n(E) — 0]. Com- 
bining the calculated n e (E) with a(E) given by Eq. ([7]) 
we obtain the required dependence of the conductivity 
a = a(n e ). 



C. Numerical results 

This subsection presents numerical results for the 
dc conductivity calculated using the time-dependent 
real space Kubo-Greenwood formalism within the tight- 
binding model. We compute the density dependence of 
the conductivity for graphene sheets with 10 and 50 lines 
in 1700 x 1000 lattice. This approximately corresponds 
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FIG. 2: (Color online) Density of states (DOS) and the rela- 
tive charge carrier concentration n e (the number of electrons 
per C atom) vs. the energy E for 10 and 50 positively-charged 
line defects described by the symmetric potential potential 
with A = 0.25u. 
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FIG. 3: (Color online) Time-dependent diffusion coefficient at different energies for 10 (a) and 50 (b) positively-charged line 
defects (A = 0.25m). 



to a relative concentration of point defects of respectively 
1% and 5%. We model the potential due to lines of 
charges by the Lorentzian function, Eq. Q , where we set 
qTFa = 0.1, which corresponds to typical electron densi- 
ties |nfP| - 5-1CT 5 atom- 1 (|raf p | - 2 ■ 10 11 cm" 2 ), see 
Fig. Q] It should be noted that gxF is no t a constant but 
weakly density dependent (<yrF oc yjn^j) . Quite remark- 
ably, the obtained results for the conductivity remain 
practically unchanged when we use different <7tf corre- 
sponding to representative electron densities considered 
in the present study, 1 • 1CP 5 < \n e \ < 5 • 10~ 5 atom -1 . 
In Appendix B we present results of more elaborated 
self-consistent calculations where we use the exact shape 
of the Thomas-Fermi potential [i.e., Eq. ^ instead of 
Eq. 0] and take into account the density dependence of 
<7tf- We found that even for a single charged line embed- 
ded in a graphene sheet, the dependence a — cr(n e ) cal- 
culated for the exact self-consistent (i.e., n e -dependcnt) 
potential ([2]) exhibits qualitatively and quantitatively the 
same sublinear behavior as in the simulations with the 
fixed Thomas-Fermi wavevector (^tfo = 0.1) and with 
Vi given by the Lorentzian function Eq. ((3|). The same 
conclusion holds for samples with 10 and 50 lines. Be- 
cause of this in what follows we will discuss the results 
for the case of the Lorentzian potential at the fixed gxF- 

Figure [2] shows the electron density n e — n(E) and 
the DOS in a graphene sheet with different number of 
charged lines. The calculated dependencies are very 
much similar to those for clean graphene and for graphene 
with a long-range Gaussian potentia h 36 ' 45 (Note that 
the DOS of graphene with short-range strong scatterers 
exhibits an impurity peak in the vicinity of neutrality 
point. ) 34 ' 36 ' 45 ' 53 ' 54 For both symmetric (not shown here) 
and asymmetric potentials, the DOS does not reach zero 
at the Dirac point, and the asymmetric potential (in con- 
trast to the symmetric one) leads to electron-hole asym- 
metry in the DOS. 

The time dependence of the diffusion coefficient at dif- 
ferent energies for the case of a symmetric potential cor- 



responding to 10 and 50 positively charged line defects 
is shown in Fig. [3] [Diffusivity curves for the case of 
asymmetric potential (not shown here) exhibit similar 
behavior.] After an initial linear increase corresponding 
to the ballistic regime, the diffusion coefficient reaches 
its maximum at t « 130 and 150 fs for 10 and 50 lines, 
respectively. These values of D = £> max are used to cal- 
culate a according to Eq. (jTj) . For times t > 150, D(t) 
decreases due to the localization effects. Similar temporal 
behavior of the diffusion coefficient was established ear- 
lier for different types of scatterers in graphene including 
long-range Gaussian and short-range potentials . 36 ' 46 ' 47 

In Fig. [4] we show the density dependence of the con- 
ductivity of graphene sheets with linear defects for the 
cases of symmetric and asymmetric potentials. The ob- 
tained dependencies show following features. 

First, the averaged conductivities exhibit a pronounced 
sublinear density dependence, see Figs. |UJc) and [4ff). 
Our numerical calculations are consistent with the recent 
experimental results for the CVD-grown graphen o 6 ' 13 ' 55 
that also exhibit sublinear density dependence. This pro- 
vides an evidence in support that the line defects repre- 
sent the dominant scattering mechanism in CVD-grown 
graphenei 6 ' 8 ' 13 Note that the calculated sublinear density 
dependence for the case of linear defects is quite different 
from the case of short- and long-range point scatterers 
where the numerical calculations show a density depen- 
dence which is close to linear j 33 i 34 i 36 ' 45 ' 56 

Second, the conductivities of samples with different 
impurity configurations exhibit significant variations be- 
tween each other, see Figs.0|a)-[3Jb), 0Jd)-|lIe). This is 
in strong contrast to the case of short- and long-range 
point scatterers where corresponding conductivities of 
samples of the same size and impurity concentrations 
practically did not show any noticeable differences for 
different impurity configurations^ 6 - We attribute this to 
the fact that in contrast to point defects, the line defects 
are characterized not only by their positions, but also by 
directions (orientations) and their intersections as well. 
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FIG. 4: (Color online) Conductivity as a function of energy E [(a),(b),(d),(e)] and relative electron density n [(c), (f)] for 
different configurations of 10 and 50 positively-charged ID defects (A = 0.25u). Conductivities in (c) and (f) are averaged over 
20 different realizations in (a), (b), and (d), (e). 
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FIG. 5: (Color online) (a) The conductivity and and relative (b) values of conductivities for two symmetrical (with respect to 
the Dirac point) energies as the functions of positive (asymmetric) potential ^~[/£ [0, A]. 



Such additional characteristics result in much more pos- 
sible distributions of the potential which, in turn, leads 
to the differences in the conductivity curves. 

Third, for the symmetric potential the conductivity 
curves are symmetric with respect to the neutrality point, 
while the asymmetric one shows the asymmetry of the 
conductivity, c.f., Figs. |U[c) and|4]T). Such asymmetry 
between the holes and electrons have been also reported 
in many transport calculations for graphene wi th p oint 
defects, for instance in Refs. 34. 42. 46. 57.581 andl3o. For 
a closer inspection of the effect of asymmetry we plotted 
the conductivities for representative energies E = ±0.5u, 



Fig. [5ja) as well as their relative differences a a iE)+a(E) 
as a function of the potential strength A, Fig. [SJb). The 
relative conductivity difference exhibits a linear behavior 
for A < 0.4u followed by saturation for larger values of 
A. A comparison of the obtained numerical results with 
the analytic predictions in the weak scattering regime 
will be given in what follows. 

We conclude this section by noting that conductiv- 
ity of large CVD-grown graphene polycrystalline samples 
with disordered grain boundaries was calculated by Tuan 
et al. 32 using the same time-dependent real-space Kubo 
method. In contrast to the long-range Thomas-Fermi 



potential considered here, the onsite potential in Ref. [32J 
is set to zero and the scattering is due to grain bound- 
aries separating domains with different crystallographic 
orientations. Even though this study did not discuss a 
functional dependence of the conductivity, a visual in- 
spection of the obtained results reveals an approximate 
linear dependence of the conductivity on the Fermi en- 
ergy, which is consistent with our results. Moreover, the 
conductance of graphene with several types of domain 
boundaries has also been shown to be a linear function 
of the Fermi energy^i We therefore speculate that the 
linear energy dependence (and thus the sublinear density 
dependence) of the conductivity is related to scattering 
off extended defects. More systematic studies of scatter- 
ing for different forms of potentials are needed in order 
to clarify this question. 



III. BOLTZMANN APPROACH 





FIG. 6: (Color online) Schematic picture showing a scatter- 
ing event and Feynman scattering diagrams considered in this 
work. The circle signifies the extended charged defect with 
charge density A and the zigzag (solid) line denotes the scat- 
tering potential (bare propagator). The transport relaxation 
rate is 0(A 2 ) in the FBA approximation (top diagram) and 
£>(A 3 ) in the SBA (bottom diagram). 



A. Formalism 

In this section we tackle the problem of dc trans- 
port in graphene with ID charged defects by means 
of semiclassical Boltzmann theory. We would like to 
stress that the full quantum calculations of Sec. [TT| and 
semiclassical kinetic theory provide complementary in- 
formation about electronic transport; while the former 
is more suitable to handle highly disordered systems or 
strong scattering regime (given practical computational 
limitations)^ semi-classical approaches yield an accu- 
rate picture of charge transport for dilute disorder and 
are often limited to the weak scattering regime (an ex- 
ception being resonant scattering which can be treated 
non-perturbatively)4i Here, the dimensionless parame- 
ter j3 = \A\L/(Kvf), with L of the order of the sys- 
tem size, defines the onset of weak scattering regime, 
i.e., /3 <C 1. Note that the simulations of the previous 
section have f3 > 10 2 , and therefore fall well inside the 
strong scattering regime. 

The effective potential of a charged line is long-ranged 
and hence we neglect intervalley scattering. Within the 
Dirac cone approximation, the semiclassical dc conduc- 
tivity of graphene at zero temperature is given by 



2h 



kFVFT{kp) . 



(8) 



In the above, the factor g = 4 accounts for spin and valley 
degeneracies, and r(kp) is the transport scattering time 
at the Fermi surface 



r(k F ) 



mv F / d9(l~ cos 9)\f(9)\ 2 



(9) 



where f(9) is the scattering amplitude at an angle 9 and 
ni stands for the (areal) density of charged lines. 

In this work we compute the scattering amplitudes in 
the second Born approximation (SBA) with respect to 



the scattering potential V(r). This allows us to improve 
over the commonly employed first Born approximation 
(FBA) by capturing the non-trivial effect of electron-hole 
asymmetry; see Fig. [5] In the Appendix C we show that 
the SBA scattering amplitude for 2D massless fermions 
is given by 



/sba(0) 




-P)x 



(u k |G (p)|u k }V(p-k)}, (10) 



where V(q) denotes the 2D Fourier transform of the 
scattering potential energy of a charged line V(q) = 
j d 2 re~ lcl ' r V(r), and Go(p) is the 2D Dirac fermion prop- 
agator for particles with energy Ep — svphkp, i.e., 



G (p) = 



1 Ep + hvF<T ■ p 



h 2 v F kp — p 



s0+ 



(11) 



The symbol s = ±1 distinguishes between electrons and 
holes, that is, s = sign(_Ep). k is the wavevector of the 
incident electron, fi,q = ft(k' — k) is the transferred mo- 
mentum (k' stands for the 'out' wavevector), 9 — Z(k',k) 
is the scattering angle, \uk) — 2~ 1 / 2 (1, se l0k ) T is the 
Dirac spinor for scattered particles, and the form factor 
E(0) = 1 + e l6 comes from graphene's sublattice sym- 
metry and precludes carriers from back-scatter. Without 
loss of generality, in what follows, we consider incident 
carriers propagating along the x-direction, k = kpe x . 

The first term inside brackets in Eq. (|10[) is propor- 
tional to the Fourier transform of the scattering poten- 
tial evaluated at the transferred momentum ftq, that is, 
the familiar FBA scattering amplitude. The remaining 
terms result from the next-order correction to the FBA 



and require the calculation of two integrals, namely 

Ij = hv F k F J ^V(k> - p)g(p)V(p - k), (12) 
X 2 = / (f^(k' - p)[»« F p • e x }g(p)V( P - k). (13) 

In writing these equations, we have defined the function 

g(p) = (Hv F )- 2 (k F -p 2 +isO + y l . (14) 

The scattering potential of an infinite line with density 
charge g = X5(x)S(z) was derived by some of the authors 
in Ref. [8| and is given by 



V(q) = 27rS(q y )- 



A 



9tf 



(15) 



where the parameter with units of energy A relates to the 
charge density of a line A according to A = ,sAe/(2e ) (in 
vacuum); note that the absolute value of A coincides with 
the definition of A as given in Sec. Ill AL The delta func- 
tion in Eq. (|15|l reflects momentum conservation along 
the direction defined by the line. For completeness, a 
derivation of this result is provided in Appendix A. 

In order to mimick the effect of lines with finite length 
we have to modify Eq. (|15|) as to allow for momentum 
transfer to occur along both spatial directions. To this 
end, we introduce a length scale associated with the line's 
average length L. In the limit of small k F L, we replace 
2ir5(q y ) — > L^- as to obtain 



V L (q) = 



LA 



m 



<7tf 



(16) 



We use this potential as a toy model for describing trans- 
port for dilute concentrations of lines of charge. The par- 
ticularly simple form of Vl (?) allows for an exact calcu- 
lation of scattering amplitudes, as shown in what follows. 



B. First Born approximation 

The FBA provides a good approximation to trans- 
port scattering rates for ID charged defects with |A| <C 
hvpL^ 1 ((3 <C 1). Within the FBA we retain only the 
first term in Eq. (1101) . The transport relaxation rate 



[TFBA(fcF)]~ 



TliV F - 



L 2 A 2 k F 



dd(l-cosd) X 



|S(6i)| 2 [Ml-cos< 



<7T F r 2 , (17) 



can be computed analytically and leads to the following 
result 



TFBA(fcp) 



1 



2k Fx /q TF (2k F + q TF ) 



vpniP 2 k F + q TF - y / q TF (2k F + q TF ) 



(18) 
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FIG. 7: (Color online) Semi-classical dc conductivity at fixed 
Thomas-Fermi wavevector as a function of Fermi energy in 
the weak scattering regime with /3 = 0.08 and qTpa = 0.01. 
Inset: The ratio Q2(oc g )/Qi(a s ) determining the amount of 
transport electron-hole asymmetry is plotted as function of 
graphene's effective structure factor a s . 



Invoking the semi-classical expression for the dc conduc- 
tivity Eq. (|2J) and the relation ^tf = 4ct g fci?, we conclude 
that a is proportional to k F , according to 



4e 2 /(4a g 
h mp 2 



-k 2 



(19) 



with f(x) = g(x)/[l + x — g(x)] and g(x) = yx{2 + x). 
The dependence of Eqs. (|18|) - (fi"9|) on the Fermi wavevec- 
tor could be anticipated from the form of the effective 
potential in Fourier space, Eq. (fTo| ; as V oc fc^ 1 , the 
relaxation time must be proportional to k F at all orders 
in the Born series, implying <j oc k F oc \n e \. In other 
words, higher order corrections to the FBA renormalize 
the mobility of carriers /j, = cr/\n e \, while preserving the 
overall dependence of a and \x on the Fermi energy. This 
property is specific to the potential Vl{q) and therefore 
is not expected to hold in models of extended charged 
defects beyond the limit of small k F L. 

We briefly discuss how the FBA conductivity compares 
with the results reported earlier in Fig. 2) The solid line 
in Fig. [7] shows the FBA conductivity at fixed Thomas- 
Fermi wavevector, i.e., qTF a = 0.01. In this case, the 
function /(4a g ) is no longer constant and the functional 
dependence of a with k F is changed to linear at small 
energies, hence resembling the Kubo results. However, 
this comparison should not be pushed too far; note that 
the slope of the FBA a versus E curve depends linearly 
on qTF, and therefore Eq. (TT9"|) (and hence the quadratic 
dependence) is recovered when the self-consistent relation 
qTF — 4:a g k F is used. In fact, 



a(k F ,q T F) 



4e 2 qTF^p 
~h niP 2 : 



k F <C ^TF 



(20) 
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whereas the Kubo simulations show a oc kp indepen- 
dently on gxF in a wide range of energies. The latter 
behavior does not occur in the weak scattering regime 
described here. 



C. Second Born approximation: electron— hole 
asymmetry 

A large sensitivity to the carriers polarity in transport 
dominated by charged lines is borne out in the numerical 
simulations of Sec. [TTJ Here, we describe this effect from 
the point of view of semiclassical transport theory. Ac- 
cording to the Fermi golden rule the transport relaxation 
rate depends on the modulus square of the scattering po- 
tential; hence, in FBA approximation, opposite charges 
±e have the same scattering amplitudes and hence can- 
not be distinguished. As observed earlier, the dependence 
of a on the carriers polarity can be captured by retain- 
ing the next term in the Born series for the scattering 
amplitude f(8) — the SB A bottom diagram in Fig. [6] 

We compute the transport electron-hole asymmetry, 
defined as 



Cdc 



'dc 



Ode 



'dc 



with a 



dc 



Cdc 



°"dc| A - 



(21) 



In the weak scat- 



tering regime, L|A| ^; hi'F, the asymmetry parameter is 
proportional to the ratio of the bottom to the top dia- 
grams in Fig. |B1 Explicitly, 



6 = 



2 signARe J d6\E(e)\ 2 (I - cos 6)11! +l 2 }V(g e ) 
J d9\E(9)mi- cos 6)\V(qeW 



where qg — 2kp sin(9/2) is the transferred momen- 
tum in elastic scattering events. Remark that 2~l(2) in 
Eq. (|2"2"1) depend on the angle through the wavevector 
k' [c.f., Eqs. CE2])-([13l)]. The derivation of this and re- 
lated results is given in Appendix D. 

Inserting the potential energy of a charged line Eq. (fTo]) 
into the above expression and performing the angular 
integration yields 



AL 



hvf 



2 (a g 



L (a g 



(23) 



The explicit form of the functions Qi(a g ) and Q2(a%) 
is given in Eqs. (15TJ) and ([B"5]) . respectively. For the 
toy model of a charged line considered here [Eq. I|16[>]. 
cross sections are proportional to kp at all orders, im- 
plying that the asymmetry parameter 5 is insensitive to 
the Fermi energy. Indeed, the electron-hole asymmetry 
depends only on the magnitude of the Thomas-Fermi 
screening through the effective graphene's structure con- 
stant, a g . In vacuum, a g ~ 2.5, and the evaluation of 
Eq. flUf yields S « 0.08 • /3. The ratio Q2/Q1 is found 
to be very sensitive to the effective screening length of 
a charged line (refer to inset of Fig. [7]); for </tf ^ kp 



(a g ^ 1) screening is very efficient and electron-hole 
asymmetry is negligible, whereas for (?tf ^ kp (a g < 1) 
the ratio Q2/Q1 can assume large values leading to an 
enhancement of the asymmetry parameter 6. 

The transport electron-hole asymmetry in scattering 
events reflects into a decrease (increase) of the SB A trans- 
port relaxation time with respect to the FBA result for 
positive (negative) Fermi energy. In fact, by expanding 
the SBA transport relaxation rate Eq. © in the small 
parameter /3, we find 



T (k F ) = [1 - sS(a g ) + (D((3 2 )] T FBA (k F ). 



(24) 



This result shows that the effect of second term in the 
Born series (bottom diagram in Fig. [5]) is to renormal- 
ize the transport relaxation time according to the carri- 
ers polarity, s, and screening strength a g . This behav- 
ior is qualitatively consistent with the numerical Kubo 
simulations (see Fig. |4l for instance). In order to make 
the comparison between the semiclassical SBA prediction 
and the simulations shown in Sec. |TT] more accurate, we 
investigate the behavior of Eq. (IM1) at fixed Thomas- 
Fermi wavevector. Note that, in this case, the asymme- 
try parameter becomes a function of the Fermi energy 
according to S — 8(q^-p/^kp). Given the behavior of the 
function Q2/Q1 at small values of its argument (see inset 
of Fig. [TJ), the asymmetry at fixed <?tf can be quite large 
even at modest kp, originating a considerable deviation 
of the conductivity at fixed (/tf from its FBA value, as 
depicted in the main panel of Fig. [JJ 



D. Comparison with Kubo simulations 

Variation of conductivity with electronic density. In 
the strong scattering regime, the numerical Kubo simula- 
tions disclose a dc conductivity that is linear in the Fermi 
wavevector, a oc kp ex In^ 1 / 2 , a very distinct behavior 
from the semiclassical prediction for the dc conductivity, 
a on kp on \n e \. At first sight, it seems that both results 
are irreconcilable; after all they focus on opposite scat- 
tering regimes. However, for the toy model of a charged 
line considered here [Eq. (|16[) ]. a oc kp at all orders in 
perturbation theory, and hence we would expect simi- 
lar semiclassical behavior even in the strong scattering 
regime. In order to investigate this question further, we 
have performed numerical Kubo simulations for a dilute 
system with a single line of charge in the strong scat- 
tering regime (see Appendix B). These simulations show 
the same functional dependence a — cr(n e ) than the sim- 
ulations of Sec. [H] for highly disordered configurations. 
This indicates a possible failure of the toy model in de- 
scribing the potential landscape of the simulations in a 
wider range of electronic densities; remark that, by con- 
struction, Eq. (|TB)) should provide a good description of 
transport only at low Fermi momentum. 

Transport electron-hole asymmetry. A decrease (in- 
crease) of the electronic mobility for electrons (holes) 
with respect to the particle-hole symmetric case V ^ 



is found in all numerical simulations with V > (Figs. 3] 
and [5]) . This effect can be ascribed to the shift of the 
charge neutrality point towards positive energy values 
caused by a potential landscape with positive sign (see 
density of states in Fig. [5J. Although the semiclassi- 
cal picture is build upon the density of states of bare 
graphene, the inclusion of higher-order diagrams (Fig. [6]) 
in the calculation of the scattering amplitude renormal- 
izes the relaxation rates according to the carriers polarity, 
thus accounting correctly for the general behavior of the 
transport electron-hole asymmetry. 



IV. CONCLUSIONS 

In this work we have considered theoretically the trans- 
port properties of graphene with extended charged de- 
fects. Recent experiments show that these defects are 
ubiquitous in chemically synthesized graphene systems 
and degrade their electronic mobilities. We modeled ex- 
tended charged defects by lines with uniform charge den- 
sities and computed their potentials according to a self- 
consistent Thomas-Fermi approach. In contrast to the 
charged point defects, the potential of a line of charge is 
screened poorly by low-energy excitations in graphene, 
resulting in long-ranged effective potentials. We con- 
sidered the regimes of weak and strong scattering by 
means of semiclassical Boltzmann theory and large-scale 
numerical evaluation of the Kubo formula, respectively. 
Whereas the semiclassical calculation reveals a famil- 
iar linear dependence of conductivity with the electronic 
density, the Kubo simulations show a robust sublinear 
dependence and conductivity nearly constant by varying 
the Thomas-Fermi wavelength by almost one order of 
magnitude. The latter is a remarkable property of ex- 
tended charged defects in graphene. 
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Appendix A: Thomas— Fermi renormalized potential 

of a charged line in graphene and respective fitting 

by a Lorentzian function 

Here we derive the effective potential of an infinite 
charged line within the Thomas-Fermi (TF) approxima- 
tion. Rearrangements of electronic density in a metal, 



around an impurity, does not alter the Fermi energy Ep, 
and thus we may write^ 



E F ~ e(r) -ecf2T>(r), 



(25) 



where e(r) and e^2D(r) are, respectively, the local energy 
of the electrons at the top of the band and the effective 
potential energy induced by the impurity charge. In our 
problem the metal is graphene (at finite densities) and 
the impurity is a charged line. Let n eq be the electronic 
density of pristine graphene, then 



e(r) = E F 



de 



dn(r) 



Mr) 



n(r)—n B 



'cqj 



(26) 



to first order in 6n(r) = n(r) — n cq . We thus arrive at 
the following relation between the potential energy and 
the charge density 



eV2D(r) 



dc 



dn(r) 
K 



2\Mr) 



n(r)=n oq 

5n(r), 



[n(r) 



(27) 



where e(r) ~ n-\/n(r) with k = UvfIcf/ y/n cq . 

The above equations show that in order to maintain 
the Fermi level constant, a change in the local electronic 
density takes place. The effective potential has to be de- 
termined self-consistently solving the Poisson's equation. 
According to the TF approximation, we have 



VVoff (r, z) = [pi, 

EqS t 



■Sp(r,z)} 



(28) 



where Eq (e r ) is & vacuum (relative) permittivity. Note 
that in the above equation tp e f((r, z ) depends on the in- 
plane coordinates r and z. We consider a line defect with 
charge per unit of length A, and orientated along the y- 
axis, 



(29) 



Pimp( 


r, z) = XS(x)S(z). 


From Eqs. (|28]>— ([29]) and 


Sp(r,z) = —e5n(r)6(z) 


2e 2 
= —y/rWf>2T>{r)5{z) , 


we arrive at the important intermediate result 


VV cff (r,z) = 

EqEi 


'2e 2 
Jn^(p 2 n (r) - \S(x) 

K 



(30) 



6(z). 

(31) 

Note that the term Sp(r) in Eq. (|2"5| is not only a self- 
consistent term, but also imposes an important geometric 
restriction by forcing the rearrangement of charge to oc- 
cur in the graphene plane. We solve the Poisson equation 
Eq. (|32j) using the Fourier transform method, viz., 



{& + <£)<PeS(.Qx,qz) 



X 

EqE t 



2qTF<f2T>(qx 



(32) 
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Effective potential, cp/U 

-- Lorentzian function, (pt/U 




FIG. 8: (Color online) The Thomas-Fermi potential (|55)l 
fitted by the Lorentzian function (137[) . Here, ^tfi = 0.1, 
U — \/(2neo£r), and the fitting parameters are A — 1.544, 
B = 0.780, C = 0.046. 



where we have defined (7tf 



/n^/ieoSrK,) = Aa g k F - 



Integrating out the q z dependence leads to 

A/(2e £r) 



P2n{qx) = / -^-p c s{qx,q z ) 



<Ztf 



(33) 



The effective potential in a real space is therefore given 
by 



<p(x) = (fi2u(x) 



X 



dq x cos (q x x) 



2£o£r Jo n 1TF 



(34) 



or, equivalently, 

ip(x) 



X 



-\ — cos (qTFx) Ci (qrFx) 

r *- 

(q T Fx) |, 



2i:e e 

+ sin (qrFx) 



2" Sl 



(35) 



where Ci and Si denote the cosine and sine integral func- 
tions. The above equation possesses the following asymp- 
totic behavior: 



A 



tp(x) 



27re:o£ : i 



\qTFX J 



A 



111 



qTFX > 1, 
qTFX -C 1. 



(36) 



The obtained expression for the effective potential, 
Eq. (|35[) . is well fitted by the Lorentzian function, 



PL (a:) 



A 



.4 



27r£ £r B + Crr 2 ' 



(37) 



where fitting parameters A, B, C can be calculated from 
the least-squares method, see Fig. [8] We use Eq. (|37| in 
the numerical calculation based on the Kubo approach. 



Appendix B: Self-consistent calculations of the 
conductivity for a single charged line 

The Thomas-Fermi wavevector ^tf entering the effec- 
tive scattering potential [Eq. @ or (|3"5"|) ] depends on the 



electron density n e . In this appendix we check how this 
dependence affects the behavior of a as compared with 
the results obtained in Sec.|II]for a fixed qtf; accounting 
for the density dependence makes our effective potential 
'self-consistent'. 

We perform our calculations as follows. In the Kubo 
method used in this study it is not possible to change the 
scattering potential while changing the energy (or den- 
sity) of the electrons. We therefore perform independent 
calculations for six different values of ^tf obtaining six 
different dependencies a = a(E) and a = a{n e ) as shown 
in Fig. IHlJa) and Fig. EJb), respectively. In each depen- 
dence a = a(n e ) we choose only one particular point (for 
both n- and p-types of charge carriers) where the electron 
density n e corresponds to (?tf used in the calculation of 
this dependence (recall that ^tf scales as gxF oc -y/|n e |). 
Combining these six points on a single plot yields a 'self- 
consistent' curve a — cr(n e ) as shown in Fig. [I3c). Fig- 
ure[9]clear demonstrates that energy and electron density 
dependencies of conductivity exhibit respectively linear 
and sublinear behaviors, which are the same as corre- 
sponding behaviors of the conductivities for the case of a 
fixed gxF, see Fig.@] Note that electron-hole asymmetry 
in Fig. IHl is weak since the source of disorder here is due 
to a single line only (c.f, with 10 and 50 lines in Fig. 2]). 



Appendix C: Scattering amplitudes in the second 
Born approximation 

The scattering problem (Hq + V — E)^-^ = 0, where 
Hq denotes the free Hamiltonian and V a potential, has 
the formal solution 



*k 



6 k + G V* k . 



(38) 



where 0k solves the free Schrodinger equation (Hq — 
E)4>k = and describes the state of the incident par- 
ticles. The resolvent is given by Gq(z) — l/(z — H ), 
where z includes an infinitesimally small imaginary part. 
In the context of the present work, Hq stands for the 
Hamiltonian of pristine graphene in the single Dirac cone 
approximation, and V refers to the potential of a charged 
ID defect (Appendix A). Although the form of V remains 
unspecified in what follows it is assumed to be a scalar 
in both sublattice and spin spaces. The spinor <pk(r) = 
(r|0k) has the form^£i^i 



k (r) = Uke' 



with 



u k 



V2 v ^ 



(39) 



(40) 



In the above, #k = arctan(fc 2/ / 'k x ) and s = siga(E). 
Switching Eq. (1381) to the position representation, we ob- 
tain the Lippmann-Schwinger equation 



tt k (r) = k (r) + f d 2 r'G (r - r')V(r')^(r') 



(41) 
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FIG. 9: (Color online) Conductivity vs. the energy (a) and electron density (b), (c) for nonself-consistent (a), (b) and self- 
consistent (c) effective potentials V ~ U — X/(2ireo£i) £ [0, A] (A = 0.25u) describing a single charged line. Curve in (c) is 
plotted combining corresponding values for a on curves in (b). The correspondence between <7tf used in calculations of the 
potential and respective n e is indicated in (a) and (b). 



where Go(r — r ') = ( r |Go(z)|r') is the Green function of 
the problem. The graphene Hamiltonian reads 



Hq = hvpcr ■ p , 



(42) 



and the Fourier transform of the Green function C?o(p) 
/ d 2 r exp (— ip ■ r) Go(r) is given by 



Go(p) = (z - hv F cr ■ p) 



(43) 



In what follows, unless stated otherwise, we set h = 1 = 
Vf ■ It is also convenient to recast Eq. (|4"3")l in the form 



G (p) = g(p)(z + cr-p) 
g(p) = (z 2 -p 2 )-\ 



(44) 
(45) 



where z = E + is0 + ; the inclusion of a small imaginary 
part is0 + amounts to consider outgoing waves (see be- 
low). For simplicity we focus on scattering of positive 
energy carriers (electrons), s = 1. We write E — k and 
evaluate the Green function in real space representation 



G (r - r') = (E- lC r-V)J -f^e^^ g(p) (46) 

(47) 



(2<r 

r(i) 



= -l(k-icr-V)H^(k\r-r'\) 

where Hn (fc|r — r'|) is the first kind Hankel function of 
order n, whose asymptotic form is that of outgoing cylin- 
drical waves. Using the property d x Hg '(x) + H{ ' (x) — 



0, the second term in Eq. (14T[) can be written in the sim- 
ple form 

a • VH^\k\r - r'|) = -kH ( ^(k\v - r'\)a e , (48) 
where we have introduced the matrix 



ere 



er u 

e e o 



(49) 



In the above, the angle 9 = 8(r, r') is defined through the 
relation (r — r') /|r — r'| = (cos9, sin6) . 



Combining Eqs. (T4"T1) - (|4"5|) we obtain the explicit form 
of the Green function of pristine graphene 



G (r-r')=-^[^ 1) (fc|r-r / |) + ^ J fT«(fc|r-r / | 
The Lippmann-Schwinger equation now reads 



(50) 



*k(r) 



ik 



rW/ 



M*)-^J d 2 r'[H^(k\v-v'\)- 

io- fl fr 1 (1) (*|r-r'|)lv(r')*k(r'). (51) 



To proceed, we assume that the main contribution to the 
scattering amplitude comes from evaluating the above 
integral within the region where |r — r'| ^> 1. We note 
that although this procedure is accurate for short-range 
potentials, yielding the exact asymptotic form of the scat- 
tered wave function, it is otherwise an approximation. 

The next step is to insert the asymptotic expressions 
for the Hankel functions 



H^(k\r-r'\)^ 



ikirlr — r f 



ik\r— r \ 



rd) 



H^(fc|r-r'|) 



D ik\r — r' 



W ikn\r — r'\ 
into the Lippmann-Schwinger equation (|5ip to get 



(52) 



(53) 



*k(r) = k (r) 



8irr 
1 + a 9 ) V(r')* k (r / ). 



''"' ' d 2 r'i" _,k ' '''' 



(54) 



In the above, we have identified the wavevector at the 
point of observation, k' = k-, and used |r— r'| ~ r— r-r'/r 
to simplify the argument of the exponentials in ([5"^j) - (l5"3"l) . 
The first term in the Born series is obtained by replac- 
ing ^k(r') —> 0k(r') = e «k in the right-hand side of 
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the Lippmann-Schwinger equation. In order to read out 
the scattering amplitude a few manipulations are still in 
order. Without loss of generality, setting 9^ — 0, and 
identifying 9^ with the scattering angle 9, we find 



= H(0)v k ', 



(55) 
(56) 



where we have defined graphene's Berry phase (form 
factor) term E(6) = 1 + e~' 10 . Substituting this result 
into Eq. (J54")) we arrive at the well-known FBA in two- 
dimensions 



with 



Mr) = Mr) + ^ 1 e^u k , 



VfU V o7T 



(57) 



(58) 



where q = k'— k is the transferred wavevector and the 
relevant constants have been restored. We note that the 
above definition yields the usual form for the scattered 
current in two-dimensions, i.e., 

J{&) = (* k (r)M* k (r)) oc l/FBA(0)|2 , (59) 

r 

with ^k(r) = ^E'k(r) — 0k (r) denoting the scattered com- 
ponent of the wave. 

We move gears to the calculation of the second term 
in the Born series. The starting point is Eq. (HIT) , which 
we iterate two times to get 

*k(r) = 0k(r) + j dV [G (r - v')V {r')fa{r')+ 



d 2 r"G (r - r')V(r')G Q (r' - r")V{v")^{r") 



(60) 



We aim to simplify the second order contribution in the 

(2) 

above expression [from now on referred to as "J^ (r)]. As 
before, we replace Go(r — r') by its asymptotic form 



It is clear that bg does not commute with the remain- 
ing terms in the integrand [remark that Gq(y' — r") con- 
tains a term proportional to a a with a = 6(r',r") ^ 6; 
c.f., Eq. (j5TJ)) ]. and hence we cannot directly identify the 
scattering amplitude as previously. Instead, we make use 
of Eq. (gBJ to write 



kk' 



/l7 9,P) / A ' V(r " ) ^ t "" 



f d 2 r'e- lk ' r 'v(r')a e (k-ia- V) e lpr ' 

or, using the definition of Fourier transform, 
d 2 p 



u k , (64) 



kk' 



(2^) : 



:V(p - k) [a e G (p)} V(k' - p)u k . (65) 



In order to identify the scattering amplitude in the second 
Born approximation (SBA) we compute the contribution 

(2) 

of ^ (r) to the scattering flux. Neglecting terms of 
fourth order in the scattering potential, we find 

Jsba(0) - (* k (r)kfl|*k(r)) = J(0) + 5J(9) , 
with J{9) given by Eq. (I5T)|) and 

^)_M)& ak , kfl | Tkk , )+ ,,. (66) 

r V o7r 



Using 

V2&0 (k + a ■ p) M k 



E(9)k + pe-*** +pe< 
E(-0)k + pe 1 ^ +pe 



i(6-<f>v) 



(157) 



where 
result 



arctan(pj,/p x ), we arrive at the following 



(uk'kelTkk') = 



d 2 p ~ 



(27T) 



V"(p-k) [jfc+pe _ ^+ 



+pe l 



ke~ lb 



V(k' 



(68) 



By the definition of scattered current Jsba(#), the SBA 
scattering amplitude is readily seen to be 



Go(r-rO^-^e*-e-*'-'^ (61) 

with bg = 1 + <jg, and insert it back into vp^ (r) as to 
obtain 



where 



r( 2 )(Vl - 



ik 



** (r) = "W Tkk '' 



(62) 



T kk ' = Jd 2 r'd 2 r"e-*' r 'agV(r')G a (r'-r")V(r")Mr") 

(63) 



f S BA(6) = J^{E(0)V(k>-k)- 



k 

8^r 
d 2 p - 



(27T) 



V{k'-p)[E(9)(k+pcos^ p ) + 



E(9)ipsincf> p g(p)]V(p-k)} 



(69) 



where we defined s(0) = s(0 + w) and dropped an in- 
nocuous phase factor — \f%. We now specialize to po- 
tentials with inversion symmetry; these potentials have 
V(q) = V"(q)* and therefore we can drop the imaginary 
term in last line of Eq. (j69|) . which is odd under the 
transformation 9 — > —9, and hence does not contribute 
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to transport cross sections. We thus arrive at our desired 
result 



/sba(0)=S(0H/— F(k'-k) 



d 2 p 

(2^) 2 



(70) 



\/(k'-p)(fc + pcos0 p ). g (p)l/(p-k)} , 
or, in a more compact form, 



/sba(0) 



5(0) 

^F^ V 87T 



t/(k' - k) 



rf 2 p 



V(k' - p)(u k \G (p)\u k )V(p - k)J ,(71) 
where ft and «f have been restored. 



Appendix D: Calculation of second Born amplitude 
for a charged line 

In this appendix we evaluate the SBA transport cross 
section for a charged line with potential given by Eq. (|16p . 
We perform an analytical calculation of the I\ contribu- 
tion [Eq. (|12j) ] and evaluate the remaining contribution 
[Eq. (|13[) ] numerically. The term X\ requires to evaluate 
the following integral 



m = 



d 2 p 



1 



(2tt) 2 k 2 
1 



i0^ 



<7tf + \Px ~ k x \ q TF + \k' x - p x \ ' 



(72) 



which we do by first performing the integration over k y 
to get 



Vi 



dp x i 



1 



1 



2tt 2p q 1F + \p x - k x \ q TF + \k' x - p x \ 



,(73) 



where po = \/k 2 — p 2 + i0 + . To proceed, we divide the 
integration range into four subintervals: p x > k x , k x > 
Px > k' x , k' x > p x > —k x and p x < —k x . Each of these 
contributions has a solution in terms of simple functions. 
We give the explicit solution for the real part of 771 . Since 
Po becomes pure imaginary for \p x \ > k, we have 




k dp x 



1 



1 1 



q-iF -Px + k x q TF + k' x - p x 
dp x 1 

k 2lT yjpl - k 2 

1 1 



9tf +Px- k x q TF - k' x + p x 

Without loss of generality we set k x = k, k' x = fccos 
The integral above then acquires the form 

1 X(0) 



(74) 



Re 771 



327ra g fcf,(l-cos0)' 



(75) 



with k ~kp and 

-2arccos(l + 4a g ) 7r + 2arcsin(l — 4a g ) 



X(0) 



ijl 



1 

2n.. 



1 



2a s 



8a g arccos (4a g — cos ff) 8a g arccosh (4a g + cos ( 



1 — (4a g — cos ( 



(4a g + cos ( 



(76) 

and where used Qtf = ^ot g kp. The remaining term to be 
computed reads 

d 2 p p x 

m 



(2?r) 2 k 2 -p 2 + i0+ 
1 1 



k x \ 



(77) 



<7tf + \K -p x \ qi F + \q x 

The explicit form of 772 is rather cumbersome and thus 
will not be given. The differential cross section is 

<r(8) = \h(8) + f 2 (9)\ 2 , (78) 

where /i(2) denotes the first (second) order contribution 
to the SBA amplitude [sec Eq. ([H])]. Defining h l( 2){9) = 
/i(2)(0)/S(0), we obtain 

a{6) = |S(6»)| 2 {h^O) 2 + 2Re [h 2 (9)} hi(6) + 0(A 4 )} , 

(79) 
and where we have used the fact that h\{0) £ R for po- 
tentials with inversion symmetry. The first term yields 
the FBA transport cross section 

r-271- 



(FBA) 
7 tp 



M(l-CQB0)\fl{0)\' 



2 „2tt 



d0 



(l-COS0)\E(d)\' 



k F f LA\ 

8tt ^ Hv F J Jo ' \2k F sin 2 (9/2) + q TF ] ' 



/LAV Qi(q) 

\ hvp J kp 



with 



Qi(a g ) 



2 * de (l-cos6)\~(9)\\ 
[2sin 2 (6V2)+4a g ] ; 



(80) 



(81) 



Remark that the transport relaxation rate is related to 
crtp according to r = (nivpat v )~ l ~ kF, and therefore 
we conclude that the dc-conductivity 



2e 2 2e 2 k F 

a = ——v F k F T[k F ) = — — — 

h h ni<j tp (k F ) 



(82) 



is a quadratic (linear) function of the Fermi wavevector 
(electronic density). The latter property is preserved at 
all orders in perturbation theory as noted in Sec. IIII1 

The second term in Eq. (|79[) yields the main correction 
to the FBA transport cross section; explicitly, 



™,2 Re[r7i(0)+v7 2 (0)] 



|H(0)|< 



2k F sm /i (9/2) + q TF 



(83) 



14 



Simplifying one obtains 



^-m 9 ^ 



with 



K) = zz 



/2l«g 



2 71 



(16 



{l-cos6)\Z{9)\ : 



(84) 



■V>sba(0), (85) 



8tt J '" 2sin 2 (6>/2)+4a g 

and V'Sba(^) = 2kpRe[r)i(6) +i]2(0)} is just a function 
of 9 and a g [recall that 171(2) varies as kp 2 ; see Eq. (|T5"|) ]. 
Finally, one obtains for the SBA transport cross section 



o, 



(SBA) 
tp 



y (^Y%ifos) + 0(A«). (86) 

^- ' \ Til 1 77- / /? T71 



n=2,3 



\flVF 



Collecting these results one obtains the following relation 
between the SBA and the FBA conductivities 



(SBA) 
(FBA) 



(FBA) 



'tp 



(LAV 



_ LA Q 2 (a s ) — 

(SBA) ^Q^r^Uw ' 

tp 



The above result shows that for A > (A < 0) the SBA 
decreases (increases) the dc conductivity with respect to 
the FBA result. Although only valid in the weak scat- 
tering regime, this dependence of the dc conductivity on 
the carrier polarity is in qualitatively agreement with the 
numerical results of Sec. [XXJ 
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